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CALCULATION  OF  THE  IONOSPHERIC  PHOTOELECTRON 
DISTRIBUTION  FUNCTION 


I.  INTRODUCTION 

Through  the  process  of  photoionization,  a fraction  of  the  solar 
radiation  which  reaches  the  earth  is  transformed  into  electron  kinetic 
energy.  Subsequent  collisions  between  the  electrons  formed  by  photo- 
ionization processes  and  the  neutral  particles  produce  the  ambient  elec- 
tron and  ion  gas  as  well  as  heat  the  neutral  atmosphere.  The  ambient 
electron,  ion,  and  neutral  gases  may  each  be  characterized  by  a Max- 
wellian distribution  at  a temperature  from  300  - ^000  °K.  Those  elec- 
trons which  are  removed  from  the  neutral  molecules  either  directly 
through  photoionization  or  through  the  collisional  formation  of  second- 
aries are  called  "photoelectrons".  They  have  energies  substantially 
greater  than  the  ones  characterizing  the  Maxwellian  electron  gas  and 
their  distribution  function  can  be  described  as  an  irregular  high  energy 
tail  on  the  Maxwellian.  The  main  source  of  heating  the  ambient  thermal 
electrons  occurs  through  their  interaction  with  the  photoelectrons. 

The  work  presented  below  describes  our  effort  to  understand  how 
energetic  electrons  are  degraded  and  transported  in  the  earth's  atmo- 
sphere. We  calculate  a distribution  function  for  the  high  energy  photo- 
electrons which  can  also  include  the  effects  of  electrons  of  even  higher 

* 

energies  which  come  into  the  atmosphere  along  magnetic  field  lines.  Our 


Note:  Manuscript  submitted  August  26,  1976. 
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goal  is  to  describe  in  detail  the  redistribution  of  the  energy  of  the 
electrons  among  the  ambient  thermal  electrons,  ions,  and  neutral  parti- 
cles. From  this  information  we  calculate  atmospheric  quantities,  such  as 
electron  temperatures  and  densities  of  excited  states,  which  are  crucial 
to  our  ability  to  understand  and  predict  the  state  of  the  ionosphere. 

Photoelectrons  can  lose  energy  through  a myriad  of  collisional 
processes  which  ionize  or  excite  the  neutrals.  For  0^  and  N^,  both  rota- 
tional and  vibrational  states  may  also  be  excited.  Excited  species  are 
important  in  the  processes  which  control  the  chemical  balance  of  the 
ionosphere  and  they  contribute  to  the  observed  dayglow  emissions.  Elec- 
tron-electron collisions  and  the  excitation  of  plasma  waves  are  respon- 
sible for  heating  the  ambient  thermal  electron  gas  and  subsequently  the 
ion  and  neutral  gas.  Any  calculation  of  the  electron  temperature  re- 
quires a description  of  the  heat  transfer  and  conduction  processes. 
Photoelectrons  are  responsible  for  the  observed  predawn  heating  of  the 
ionosphere  due  to  the  influx  of  electrons  along  magnetic  field  lines. 

This  occurs  when  the  conjugate  ionospheric  region  is  in  daylight.  If  we 
wish  to  properly  account  for  all  of  these  processes,  we  need  to  know  the 
details  of  the  atomic  and  molecular  structure,  the  relevant  interaction 
cross  sections,  and  the  plasma  interactions. 

Photoelectron  transport  calculations  have  been  made  by  Nisbet 
(I968)  and  Swartz  (1972),  Nagy  and  Banks  (1970)  and  Banks  and  Nagy  (1970), 
Cicerone  and  Bowhill  (I970,  I97I),  and  most  recently  by  Mantas  (1975). 

The  Nisbet  method  is  based  on  diffusion  theory  and  determines  the  photo- 
electron flux  by  solving  a set  of  coupled  flux  equations.  Escaping 
electrons  are  accounted  for  by  including  an  additional  loss  term  which 


is  separately  calculated.  Swartz  modified  this  approach  by  directly 
including  the  flux  divergence  term.  Nagy  and  Banks  solved  for  the  elec- 
tron flux  in  a two-stream  approximation  and  thus  attempted  to  treat  the 
high  altitude  transport  properties  of  photoelectrons.  Cicerone  and 
Bowhill  used  a Monte  Carlo  method  which  follows  the  path  of  the  individ- 
ual particles.  The  limitations  and  problems  inherent  in  each  of  these 
three  approaches  have  been  analyzed  in  detail  by  Cicerone  et  al.  (1973). 
The  Nisbet  and  Swartz  method  requires  the  calculation  of  a diffusion 
coefficient  for  electrons  through  the  atmosphere  which  depends  on  the 
pitch  angle  distribution  and  on  the  details  of  the  collisional  scatter- 
ing. A two-stream  model  such  as  Nagy  and  Banks'  can  omit  the  essential 
features  of  angular  anisotropies  in  the  particle  distributions.  A Monte 
Carlo  method  is  in  principle  the  most  comprehensive  and  does  not  require 
any  information  a priori  on  the  pitch  angle  distribution.  It  is  limited, 
however,  by  its  computer  cost  and  core  requirements. 

The  transport  equation  approach,  such  as  that  used  by  Mantas 
(1975)  photoelectrons  and  by  Strickland  et  al.  (1976)  for  auroral 
electrons,  is  currently  the  most  feasible  and  versatile  method.  Below 
we  show  solutions  of  the  Boltzmann  equation  which  give  the  photoelectron 
distribution  function.  As  a biproduct  we  find  the  population  density 
of  excited  ionic  and  neutral  states  from  which  airglow  emissions  can  be 
calculated.  With  minimal  additional  effort  the  effects  of  high  energy 
particle  fluxes,  electric  and  magnetic  field  effects,  and  varying  solar 
conditions  can  be  included  in  the  calculation.  The  numerical  method 
is  based  on  the  one  described  by  Strickland  et  al.  (I976)  for  auroral 
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electrons  and  will  be  discussed  in  more  detail  in  a later  section.  In 
this  paper  we  will  describe  model  calculations  of  the  photoelectron  dis- 
tribution function  and  give  some  comparisons  with  observed  spectra. 
Specific  studies  using  these  results,  such  as  electron  temperature, 
airglow  calculations,  or  conjugate  point  effects  are  reserved  for  later 
papers. 


II.  THEORETICAL  BACKGROUND 

The  photoelectron  distribution  function,  f(r,  v,  t),  is  found  from 
the  solution  of  the  Boltzmann  equation, 


— + v-V  f + v.V  f = Q(r,  v,  t)  + — , (1) 

at  " r " v “ at 

c 


where  df/at  represents  the  local  variation  of  f with  time,  the  second 
term  describes  the  change  in  f due  to  spatial  gradients,  the  third  term 
describes  the  presence  of  external  forces,  Q is  the  local  production 
rate  of  primary  electrons  due  to  solar  radiation,  and  af/^t|  reflects 
the  change  in  f due  to  elastic  and  inelastic  collisions  with  other  elec- 
trons, ions,  and  neutral  constituents.  For  our  present  purposes  we  have 
introduced  a number  of  simplifying  assumptions  which  are  described  in 
the  paragraphs  below. 

First,  we  assume  that  photoelectrons  are  transported  predominately 
along  magnetic  field  lines  and  that  the  atmosphere  is  horizontally 


stratified  along  these  field  lines, 
at  middle  and  high  latitudes,  allows 


This  approximation,  which  is  good 
us  to  assume  azimuthal  symmetry 


k 
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about  a magnetic  field  line  and  also  to  write  the  spatial  gradient  in 


as  a derivative  along  a direction  z which  is  parallel  to  B.  We 


then  need  to  consider  transport  only  along  z.  The  time  derivative  of 


the  velocity,  v,  is  proportional  to  the  electric  field,  E,  and  to  v x B 


where  e is  the  charge  on  an  electron  and  m is  its  mass.  Since  we  con- 


sider transport  only  along  z and  assume  no  electric  fields  are  present 


this  term  does  not  contribute  to  Eq.  (1) 


Since  interaction  with  the  photoelectrons  is  the  primary  source 


of  heating  the  ambient  thermal  electrons , we  now  consider  how  to  include 
this  in  Eq.  (l).  Butler  and  Buckingham  (1962),  Itikawa  and  Aono  (1966), 
and  Perkins  (1965)  have  looked  in  some  detail  at  the  way  in  which  a 


charged  particle  loses  energy  as  it  moves  through  a thermal  plasma  When 


we  apply  their  work  to  the  ionosphere,  we  conclude  that  energy  loss  to 


ions  is  negligible  and  that  we  only  have  to  consider  photoelectron- 


thermal  electron  interactions.  Energy  loss  is  due  to  the  combined 


effects  of  two-body  Coulomb  collisions  and  the  collective  effect  of 
Cerenkov  emission  of  plasma  waves.  Schunk  and  Hays  (1971)  have  expressed 


the  energy  loss  of  a high  energy  electron  in  a lower  energy  electron  gas 
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where  E is  the  photoelectron  energy, 


v = 


v; 


2 kT 


(b) 


is  the  photoelectron  velocity,  cu^  is  the  plasma  frequency 


= 


m 


(5) 


n^  and  T are  the  thermal  electron  number  density  and  temperature,  re- 
spectively, and  -Ln  V is  Euler's  constant.  The  second  formula  in  Eq.  (3) 
describes  energy  loss  for  those  energies  in  which  the  deBroglie  wave- 
length is  greater  than  the  classical  distance  of  closest  approach  and  a 
quantum  mechanical  treatment  is  required.  Equation  (3)  is  incorporated 
into  Eq.  (1)  by  the  method  described  by  Mantas  (1975)>  who  treats  this 
loss  process  as  an  extra  "friction"  term  which  is  added  to  v in  Eq.  (2). 
Note  that  we  are  assuming  a continuous  energy  loss  of  the  photoelectrons 
to  the  ambient  electrons  and  that  the  photoelectrons  are  not  deflected 
in  this  process.  We  have  included  this  process  In  the  third  term  in 
Eq.  (1). 
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We  wish  to  solve  for  the  flux,  0(t,  z,  E,  u),  which  is  related  to 
the  distribution  function  through 

0(t,  z,  E,  u)  = v(e) f ( t , z,  E,  u)  . (6) 

Then  for  the  assumptions  we  have  just  discussed,  Eq.  (1)  becomes 

- — + U — + Un  {(e)[  — — \ = ••*,  (7) 

v 3t  ?z  e V 2 E ?)E  / 

where  u,  is  the  cosine  of  the  pitch  angle,  the  angle  between  v and  B. 

The  function  ,£(e)  is 


where  E « l4  eV. 

o 

In  the  work  described  below  we  assume  that  the  steady  state  dis- 
tribution function  gives  an  adequate  description  of  the  photoelectron 
flux  and  set  ?»0/?)t  = 0.  This  is  justified  by  the  short  decay  time  of 
the  initial  distribution  compared  to  any  changes  in  the  source  term 
which  produces  primary  electrons.  Such  an  assumption  is  even  reasonably 
accurate  over  the  sunset  period,  except  in  a small  altitude  range  from 
which  the  sun  is  disappearing  at  that  time.  We  also  transform  Eq.  (7) 
to  solve  for  0 as  a function  of  t,  the  electron  scattering  depth,  which 


1 
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where  n^z)  is  the  number  density  and  cr  (e)  is  the  total  scattering  cross 
section  for  species  l.  With  these  changes,  Eq.  (7)  becomes 


where  U and  U are  the  cosines  of  the  pitch  angles  associated  with  the 


incident  and  outgoing  directions  and  E and  E are  incident  and  outgoing 


electron  energies.  The  function  x gives  the  probability  of  scattering 


from  the  initial  to  the  final  energy  and  direction.  It  can  be  written 


JTSSS 


The  collision  term  in  Eq.  (1),  5)f/bt|c>  appears  now  as  the  first 
and  third  term  on  the  right  hand  side  of  Eq.  (10).  The  first  term  is 
proportional  to  the  number  of  particles  leaving  a volume  element  in  phase 
space  and  the  third  term  describes  how  particles  are  scattered  into  a 
volume  element.  We  consider  collisions  between  photoelectrons  and 
neutral  particles  which  we  assume  are  in  thermal  equilibrium.  Collisions 
with  ions  cause  negligible  loss  of  energy  and  electron  collisions  have 
already  been  included.  The  redistribution  function,  R,  is  related  to 
the  cross  sections  by 

D U,  E',  E) 

R}(u',  fi,  E',  E)  = -J z , (13) 

1 ct  ( E) 

t 

where  the  sum  is  over  all  processes  and  a (E)  is  a total  cross  section 
for  the  species  t.  The  scattering  processes  which  we  include  are 
elastic  scattering  by  neutrals  and  excitation  and  ionization  of  ground 
state  neutral  particles.  We  consider  three  neutral  species,  Na,  02,  and 
0 and  the  relevant  cross  sections  are  summarized  in  Appendix  A. 

The  source  function,  Q(T,  E,  y) , is  the  number  of  electrons  pro- 
duced directly  by  solar  photoionization.  The  important  processes  con- 
tributing to  Q are  included  and  the  photon  fluxes  used  are  those  given 
by  Hinteregger  (1970)  and  those  tabulated  by  Donnelly  and  Pope  (1973). 

The  relevant  portion  of  the  solar  spectrum  has  been  divided  into  90 
bands  and  lines.  More  calculational  details  and  the  relevant  ionization 
cross  sections  are  summarized  in  Appendix  B. 
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where  the  indices  i and  n refer  to  pitch  angle  and  energy,  respectively, 
j runs  over  pitch  angles,  and  t runs  over  species.  We  have  assumed  that 
the  flux  is  linear  in  fj,  and  quadratic  in  In  E over  the  E,  n grid.  At 
the  highest  energy  level,  El5  the  flux  has  a power  law  dependence 

0(T,  E > E1}  |u)  = 0(T>  Elf  4)  y > ( 1 5) 

where  a is  chosen  to  provide  a smooth  continuation  of  the  incident  flux 
above  Ex.  The  source  function,  S,  includes  those  primary  electrons 
formed  directly  by  photoionization  which  are  described  by  Q as  well  as 
those  created  through  the  degradation  of  high  energy  electrons.  A 
detailed  description  of  this  matrix  representation  is  given  in  Strickland 
et  al.  (I976). 

The  integration  routine  CHEMEQ  by  Young  and  Boris  (1973)  bas  been 
used  to  integrate  Eq.  (14) . For  a given  energy,  we  begin  by  specifying 
the  downward  flux  at  r = 0 and  setting  the  upward  fluxes  equal  to  those 

10 


produced  directly  by  solar  photoionization  in  the  upward  direction. 
Integration  proceeds  for  downward  fluxes  from  t - 0 to  t , the  depth  at 
the  lower  boundary.  The  solutions  are  stored  over  a preselected  T-grid. 
Starting  with  an  upward  flux  of  zero  at  tq,  integration  is  then  performed 
for  upward  fluxes  from  tq  to  0 using  as  downward  fluxes  the  previously 
stored  values.  The  upward  fluxes  are  stored  on  the  same  T-grid  used  for 
storing  downward  values.  This  completes  the  first  of  several  iterations 
needed  to  achieve  convergence.  Through  a technique  which  projects  the 
iterated  solutions  to  their  final  values,  as  few  as  three  or  four  itera- 
tions achieve  excellent  convergence  to  the  final  au:wer.  The  source 
term  Q is  evaluated  by  using  a solar  deposition  routine  written  by  T.  R. 
Young  and  described  in  Oran  et  al.  (197^). 

The  FORTRAN  routines  which  solve  Eq.  (l4)  have  been  run  on  a 
360/91  computer.  For  a test  case  of  50  energies  (.5  to  120  eV)  and  10 
pitch  angles,  we  require  about  2 minutes  to  evaluate  all  of  the  R-mai rices 
for  the  three  neutral  species.  These  are  stored  since  they  need  to  be 
evaluated  only  once  for  a specified  E and  u grid.  The  final  program, 
PHOTO,  uses  the  R-matrices  and  the  primary  electron  data  to  solve  for 
the  photoelectron  fluxes.  This  requires  another  two  minutes  or  less  of 
computer  time  and  requires  ~ 4-50  K-bytes  of  storage.  The  times  given 
include  those  required  for  a sizable  amount  of  I/O.  The  program  is 
flexible  enough  so  that  both  AE  and  AU  can  vary  throughout  the  range 
of  energies  and  pitch  angles.  Thus  we  can  easily  look  at  one  particular 
energy  range  on  a much  finer  grid  scale.  Furthermore,  if  there  are  core 
storage  problems  the  program  can  be  run  in  tandeir^  where  calculations 
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from  higher  energies  are  used  as  initial  conditions  for  lower  energy 
calculations . 

In  the  next  section  we  define  a "standard"  calculation  for  speci- 
fied physical  conditions  and  show  how  it  is  affected  by  changing  physical 
parameters  as  well  as  parameters  peculiar  to  the  numerical  scheme  and 
choice  of  grids.  As  an  initial  test  of  the  method,  we  compare  a typical 
calculation  with  rocket  and  satellite  measurements  of  photoelectron 
fluxes.  The  calculations  in  the  next  section  look  in  detail  at  fluxes 
below  100  eV.  We  note,  however,  that  the  numerical  scheme  is  capable 
of  treating  much  higher  energy  electrons,  such  as  those  in  the  keV  range 
describing  auroral  electrons  (Strickland  et  al.,  197o) . 

IV . RESULTS 

We  have  two  goals  in  presenting  the  calculations  in  this  section. 
First,  we  wish  to  show  that  the  theoretical  and  numerical  techniques 
outlined  in  Sections  II  and  III  produce  realistic  results.  We  do  this 
by  choosing  a "standard"  calculation  and  showing  that  it  is  in 
good  agreement  with  observed  photoelectron  fluxes.  Second,  we  show  the 
sensitivity  of  our  calculations  to  variations  in  the  important  physical 
parameters  and  also  to  the  internal  parameters  controlling  the  numerical 
calculation.  More  detailed  comparisons  to  observations  and  further 
studies  using  these  results  will  be  described  in  subsequent  papers. 

A.  Standard  Calculation  --  Qualitative  Comparison  to  Measurements 
Our  "standard"  calculation  is  for  a zenith  angle,  x,  equal 
to  48°  and  a neutral  atmosphere  which  is  described  by  a Jacchia  (1^71) 

model  with  an  exospheric  temperature,  T , equal  to  Q00  °K,  The  neutral 

00 
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number  densities  and  the  assumed  ambient  electron  density  are  shown  in 
Fig.  1.  The  primary  electron  spectra,  which  result  directly  from 
photoionization  by  solar  radiation,  is  typical  of  low-medium  solar  con- 
ditions. We  have  shown  this  in  Fig.  2 for  several  altitudes.  For  this 
photoelectron  calculation,  we  have  assumed  that  the  primary  flux  is  pro- 
duced isotropically  in  pitch  angle.  In  the  standard  calculation  we 
have  evaluated  the  equilibrium  photoelectron  fluxes  at  2.5  eV  intervals 
from  100  eV  to  10  eV  and  at  2 eV  intervals  below  this  energy.  At  each 
energy  and  altitude,  the  flux  is  calculated  for  ten  evenly  spaced  values 
of  the  cosine  of  the  pitch  angle  from  1 (down  the  field  lines)  to  -1 
(up  the  field  lines).  The  values  of  u do  not  have  to  be  evenly  spaced 
and  can  vary  depending  on  the  resolution  we  require.  The  spacing  in 
altitude  is  exponential  and  we  have  arbitrarily  chosen  to  look  at  115 
to  750  km.  For  all  of  the  calculations  presented  below,  the  incoming 
flux  at  the  highest  altitude  has  been  set  to  zero.  The  effects  of  high 
energy  particle  precipitation  will  be  discussed  in  another  paper. 

Figures  5a- f compare  measured  and  the  standard  calculated 
photoelectron  fluxes  for  selected  altitudes  from  ~ 120  km  to  ~ 600  km. 

The  vertical  axes  show  the  total  flux  divided  by  4rr.  No  attempt  has 
been  made  here  to  reproduce  the  neutral  atmosphere,  electron  density,  or 
the  solar  conditions  characteristic  of  any  particular  dataset.  The 
observations  are  those  of  Knudsen  and  Sharp  (1972),  Galperin  et  al.  (1972), 
and  Doering  et  al.  (1970,  1975,  and  1976).  We  will  not  give  a detailed 
evaluation  of  the  accuracy  of  the  data,  though  we  note  that  the  1970  and 
1972  observations  are  less  accurate  at  low  altitudes  and  low  energies 
due  to  the  photon-produced  background  counts.  The  later  Doering  et  al. 


fluxes  (1975,  1976)  are  the  most  accurate  measurements  although  they  are 
inaccurate  at  energies  above  60  eV  due  to  stray  photons  and  cosmic  ray 


counts.  The  dip  in  observed  fluxes  in  the  low  energy  range  (Doering 


1975)  is  shown  in  Fig.  j>d.  Doering  et  al . (1976)  has  explained  that 


this  is  erroneous  and  the  most  recent  data,  such  as  shown  in  Fig.  3c 


for  180  - 190  km,  agrees  much  better  with  the  theoretical  calculations 


We  note  that  the  shape  and  magnitude  of  the  flux  at  ~ 250  km  agrees 


quite  well  with  Victor's  calculations  (1975)  which  do  not  include 


Below  ~ 250  km, the  photoelectron  flux  is  isotropic  in  pitch 


angle  for  all  energies.  At  these  low  altitudes  collision  frequencies 


are  high  enough  so  that  any  initial  directional  preference  is  soon 


eliminated.  Figures  4a-d  show  the  flux  as  a function  of  pitch  angle  for 


selected  energies  at  a fixed  altitudes.  At 


moving  up  the  field  line  can  be  several  orders  of  magnitude  larger  than 


the  downward  flux.  It  is  at  these  F- region  altitudes  and  above  where 


both  transport  and  anisotropy  effects  are  important  in  a description  of 


photoelectrons.  Figures  5a  and  b show  the  variation  of  the  flux  as  a 


function  of  pitch  angle  at  several  altitudes  for  20  and  40  eV  electrons 


For  completeness,  we  show  total  flux  vs  altitude  for  fixed  energies  in 
Fig.  6 and  total  upward  and  total  downward  flux  vs  altitude  for  fixed 


B.  Numerical  Sensitivity  Tests 


The  calculations  described  below  show  the  effect  of  varying 


parameters  which  describe  the  physical  environment  as  well  as  those  which 


control  the  details  of  the  numerical  method. 

The  effect  on  the  photoelectron  flux  due  to  an  increase  in 
zenith  angle  is  most  important  at  low  altitudes.  Figures  8a  and  b com- 
pare the  calculated  flux  for  three  zenith  angles,  X = 0°,  48°,  and  90° 
at  119  and  5^5  km.  The  X = 0°  and  the  X = 90°  calculations  have  assumed 
Jacchia  (1971)  neutral  atmospheres  characterized  by  T = 1100°  K and 
870c  K respectively.  The  y = 48°  neutral  atmosphere  was  shown  in  Fig.  1 
and  it  corresponds  to  T = 900°K. 

CO 

The  thermal  electron  gas  heating  rate  is  found  directly  from 
the  photoelectron  flux  using 


Vz) 


/max 

dE 

ds 


$ (E,z)dE  , 


where  dE/ds  is  given  by  Eq.  (3)  and  $ as  in  the  total  photoelectron 
flux  at  altitude  z.  The  basic  requirement  on  Eq,  the  lower  limit  of  the 
integral,  is  that  Eq  » kT,  where  T is  the  ambient  electron  temperature. 
In  the  present  calculation  we  have  set  Eq  = 1 .0  eV.  This  heating  rate 
is  shown  in  Fig.  9 for  the  three  zenith  angles  mentioned  above. 

The  assumption  of  isotropic  elastic  scattering  is  common  in 
photoelectron  calculations  and  we  have  made  several  tests  to  evaluate 
its  validity.  The  standard  calculation  assumes  that  the  scattering  of 
primary  electrons  is  forward  peaked  (see  Strickland  et  al.,  I976,  for 
details) . An  isotropic  distribution  is  assigned  to  the  secondaries 
produced.  For  altitudes  below  250  km,  any  anisotropies  in  the  flux  are 
removed  by  the  high  collision  frequencies.  Above  this  altitude,  the 
effects  of  assuming  purely  isotropic  scattering  are  apparent  in  the 


magnitude  of  the  total  flux,  but  more  pronounced  in  the  pitch  angle 


distribution.  Figures  ICaand  b,  and  Fig.  11a  and  b compare  the  standard 
calculation  with  one  that  assumes  isotropy  in  elastic  scattering.  As 
we  would  expect,  the  photoelectron  fluxes  show  less  variation  in  pitch 
angle  when  isotropic  scattering  is  assumed. 

We  now  consider  the  resolution  in  pitch  angle  which  is  required 
to  get  accurate  solutions.  We  have  calculated  photoelectron  fluxes  with 
three  different  grids  that  contain  k,  10,  and  20  points.  For  an 
energy  grid,  we  have  used  a 20  energy  point  grid  with  iE  = 5 eV  from 

90  to  5 eV.  The  results  are  shown  in  Figs.  12  and  13,  We  conclude  that 

changing  the  fj,  grid  will  not  substantially  alter  either  the  shape  or 
magnitude  of  the  averaged  flux.  From  Figs.  13  we  see  that  the  *4-  point 

grid  (4  = + 1,  0.5,  - 0.5,  - 1)  gives  essentially  the  same  answers  at 

these  values  of  as  the  10  or  20  point  grids,  although  the  points  are 
too  far  apart  to  construct  the  true  shape  of  flux  vs  when  anisotropies 
occur.  The  10-point  grid  does  give  essentially  identical  answers  to  the 
finer  grid. 

The  standard  calculation,  which  has  a finer  energy  grid 
(^E  = 2.5  eV)  and  10  points  in  the  jx  grid,  is  also  shown  on  Figs.  12  and 
13.  At  low  altitudes  the  magnitude  of  the  curves  will  be  altered 
significantly  by  greater  energy  resolution.  Figure  12b  shows  more 
structure  appearing  as  the  energy  resolution  is  increased.  The 
differences  in  magnitude  above  60  eV  are  due  to  the  different  limits 
the  calculation.  The  standard  was  calculated  to  120  eV  to  insure  accuracy 
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The  20  energy  test  calculations  were  only  to  90  eV.  A 100 


point  energy  test  differs  from  the  50  energy  standard  calculation  by  a 


indicating  that  the  50  energy  calculation  is  accurate 


maximum 


enough  to  give  the  correct  magnitude  and  energy  for  the  flux 


V.  CONCLUSION 


We  have  presented  the  details  of  our  calculation  of  the  iono 


spheric  photoelectron  distribution  function.  Since  we  solve  the 


Boltzmann  equation,  non-local  effects  are  included  and  the  results  are 


valid  throughout  the  D-,  E-,  and  F-  regions.  The  method  of  solution  is 


computationally  fast  and  accurate  so  that  a number  of  important  iono 


spheric  studies  are  possible.  Future  work  will  look  at  phenomena  such 


as  high  altitude  escape  fluxes,  airglow,  electron  temperatures,  and  the 


effects  of  incoming  electron  fluxes  in  the  ionosphere 
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Appendix  A 


Electron  Impact  Cross  Sections 


The  elastic  cross  sections  for  N , 0 and  0 are  shown  on  Figs. 

2 2 

14,  15,  and  16,  respectively.  For  0,  the  total  and  differential  elastic 
scattering  has  been  taken  from  Blaha  and  Davis  (1975).  The  values  of 

Shyn  et  al.  (1972)  have  been  used  for  N . The  0 values  have  been  obtained 

2 2 

from  Trajmar  et  al.  (1972),  Watson  ( 1967)  and  Dehmel  ( 197o) . These 
cross  sections  are  consistent  with  those  found  by  subtracting  the  total 
inelastic  cross  section  from  the  measured  total  scattering  cross  section. 
The  method  by  which  the  differential  inelastic  cross  section  is  in- 
corporated in  the  model  has  been  discussed  in  Strickland  et  al,  ( 197°) • 

References  are  given  in  Tables  I,  II,  and  III. 


Table  I - Inelastic  Cross  Sections  for  N 


State 

Reference 

3E 

Winters  (l 966);  P.  S.  Julienne  (private  communication) 

a 3i+ 

u 

Borst  (1972) 

b 3ng 

Stanton  and  St.  Johns( 19^9) 5 Chung  and  Lin  (1972) 

c3nu 

Finn  et  al.  (1972);  Aarts  and  deHeer  (1969) 

Vibration 

Schulz  (1964);  Engelhardt  et  al.  (1964) 

w3au 

Chung  and  Lin  (1972);  normalized  to  Chutjin  et  al. 
(1975) 

a xn 

g 

Green  and  Stolarski  (1972) 

b X 

Green  and  Stolarski  (1972) 

b aE  + 

u 

Green  and  Stolarski  (1972) 

Rydbergs 

Green  and  Stolarski  ( 1972) 

Ionization 

. --  - - 

Rapp  and  Englander-Golden  ( 1965) 

IPS* 
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Table  II  - Inelastic  Cross  Sections  for  0 


Reference 


Schulz  and  Dowell  (19^2);  Hake  and  Phelps  (1967); 
Trajmar  et  al.  (1971);  Trajmar  et  al.  (1972); 
Wong  et  al.  (1973) 

Trajmar  et  al.  (1971) 


b XE  + 


Trajmar  et  al.  (1971) 


A 3I  + 


Watson  et  al.  (1967) 


Trajmar  et  al.  (1972) 


"9.9"  Feature  Watson  et  al.  (1967) 


Rydbergs 


Watson  et  al.  (1967) 


Ionization  Watson  et  al.  (1967);  Silverman  and  Lassettre  (1957); 
Kieffer  and  Dunn  ( 1966)  Rapp  et  al. 


21 


Table  III  - Inelastic  Cross  Sections  for  0 


State 

Reference 

1D 

Henry  et  al.  (1969) 

XS 

Henry  et  al.  (1969) 

5S 

Davis  et  al.  (1975);  Julienne  and  Davis  ( 197^) 

3S 

Davis  et  al.  (1975);  Julienne  and  Davis  ( 197&) 

(5s)  3d 

Davis  and  Blaha  (1976) 

Rydberg 

(3s)  XD 

Davis  and  Blaha  ( 197^) 

Calculation  of  Primary  Electron  Spectrum 


The  electrons  produced  directly  by  solar  photoionization  of  the 


ambient  neutral  particles  in  the  atmosphere  will  be  referred  to  as 


expression 


ion 

■ 

i,  i 


gives  the  number  of  primaries  produced  at  altitude  z by  solar  radiation 


energy  of  the  primaries  produced  is  E - I.  .,  where  I.  is  the  ionization 

Xj  y 1 Xj  y 1 

potential  of  ionization  state  i of  neutral  species  i,  cr«°?  (e)  is  the 


corresponding  ionization  cross  section,  and  n.(z)  is  the  density  of  f, 


The  solar  radiation,  characterized  by  the  flux 


is  attenuated  as  it  penetrates  the  earth's  atmosphere.  This 


total  absorption  cross  section  of  neutral  species  m and  the  integral  is 


taken  from  the  altitude  point  z along  the  line  of  sight  to  the  sun 


The  atmospheric  species  are  ionized  by  radiation  with  energies 
down  to  1026  i..  Characteristic  fluxes  for  these  energies  for 


various  solar  conditions  have  been  culled  from  Hinteregger  (1970) 


Donnelly  and  Pope  (1973)>  Swider  (1969),  Banks  and  Kocharts  (1975)  and 
Smith  and  Gottlieb  (197^0-  Total  absorption,  ionization,  and 


dissociative  ionization  cross  sections  in  the  energy  range  1 - 1026  A 


have  been  tabulated  from  the  references  given  in  Table  IV.  We  have 


included  ionization  from  the  neutral  ground  state  to  a number  of  singly 


positive  excited  states  of  the  ion.  These  have  been  listed  in  Table  V. 


A typical  primary  spectra  is  shown  in  Fig.  2, 


mm mwaum  * 


Table  IV 

References  for  Photoionization  and  Absorption  Cross  Sections 


Species 

References 

Comments 

N 

2 

Hubbell  (1971'' 

X-rays:  cr(Nn)  « 2 x 0 (N) . 

Lee  et  al . ( 1975) 

Total  absorption,  l80  - 700  Z. 

Cook  and  Metzger  (1964) 

Total  absorption  and 

ionization,  oOO  - 1000  l . 

Wight  et  al.  ( 197°) 

Dissociative  ionization 
rates,  20  - 60  eV. 

Stolarski  and  Johnson  (1972) 

General  compilation 

Hudson  (1971) 

Compilation  and  evaluation. 

0 

2 

Hubbell  (1971) 

Lee  et  al.  (1973) 

X-rays:  c(02)  *2  xa(0). 

Blake  and  Carver  ( 19 '5 5 ) 

Partial  photoionization 
cross  sections. 

Cook  and  Metzger  (1964) 

Total  absorption  and 

ionization,  6000  - 1000  A. 

Matsunaga  and  Watanabe  ( 1967) 

Absorption  and  photoionization, 

580  - 1070  A . 

Fryar  and  Browning  (1973) 

Dissociative  ionization 
for  304  and  584  A. 

Stolarski  and  Johnson  (1972) 

General  compilation 

Hudson  (1971) 

Compilation  and  evaluation. 

0 

Hubbell  (1971) 

X-rays . 

Henry  ( 1967,  1970) 

Partial  photoionization 
cross  sections. 

mamrw- 
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- Neutral  atmosphere  and  electron  density  model 
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Fig.  2 ~ Primary  electron  production  rate  at  selected 
altitudes  for  standard  calculation 
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Fig.  3 — Comparison  of  the  photoelectron 
to  measured  fluxes  at  i 
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Photoelectron  fluxes  as  a function  of  u,  the  cosine  of  the  pitch  angle, 
for  selected  altitudes  at  fixed  energies 
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. 8 — Comparison  of  photoelect*-on  fluxes  as  a function  of  energy 
for  x = 0°,  48°,  and  90°  at  selected  altitudes 


lermal  electron  heating  rates  as  a function  o: 
altitude  for  X = 0°,  48°,  and  90° 
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Fig.  10  — Photoelectron  flux  as  a function  of  energy  for  isotropic 
and  forward  peaked  elastic  scattering 
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st  of  angular  resolution.  Photoelectron  flux  as  a function  o: 
the  cosine  of  the  pitch  angle,  at  selected  energies. 
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Fig.  15  — Electron  impact  excitation  cross  sections  for  02 
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